Post detection selection combining diversity receivers for mobile and indoor radio channels

ABSTRACT

A method of estimating the quality of a communication channel from a differential phase angle between a received signal and the corresponding transmitted phase angle employs determining a signal to impairment ratio (SIR) as an indicator of channel state information (CSI). A maximum likelihood estimation procedure is employed to calculate this CSI metric as a function of the differential phase angle between be received signal and transmitted signal An alternate embodiment employs an estimation that incorporates average SIR information for a Rayleigh fading channel. Since CSI is derived from the phase angle of the received signal, and does not require signal amplitude information, it is attractive for use with differential detectors, phase-locked loops (PLLs) and hard-limited signals. The CSI provided can be used for implementing post detection selection diversity, by selecting the signal from a plurality of antennae which has the best SIR.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention related generally to communication of digital information and more specifically to channel state information estimation employed in decoding transmitted digital information.

2. Description of Related Art

Digital information is transmitted from one place to another through a channel. A channel is a general term referring to the medium through which a signal is transmitted. The channel may be a radio transmission, local area network transmission, telephone transmission, etc. In each case, the received signal differs from the transmitted signal due to the effects of the transmission through the channel. In most digital communication systems it is desirable to have knowledge of this channel information, known as channel state information (CSI) at any instant of time in order to properly decode the transmitted information.

Typically, CSI is derived using synchronization ("sync") signals, which may be sounding signals or a set of pilot symbols sent by the transmitter. The sync signals are predetermined and stored in the receiver. During transmission, the receiver compares the stored sync signals with the received signals and estimates the channel impulse response (CIR). This may be performed by a mean squared error or a least squares error formulation.

The received signal may be viewed as a vector comprising an amplitude component and a phase angle. In many applications only the phase angle of the received signal is retained for decoding the signal, as in the case of differential detection. The motivation for retaining only the phase angles is that the receiver may be simplified. Also, with differential encoding and detection, the transmission is more resistant to fading and Doppler frequency shifts. In other situations the received signal may be hard-limited, where the amplitude is limited to a maximum amplitude value, and hence only the phase angles of the received signal are available for decoding. The transmitted signals is sent with message data interleaved in periods between the sync signals. Therefore, in the absence of known sync signals (pilot or sounding), the CSI must be estimated from previous CSI or determined from the received signal directly.

CSI may be used in a number of ways to enhance detection of the transmitted signal. Some examples are:

Post-Detection Selection Diversity For Radio Transmission

A typical implementation of post-detection selection diversity will employ several antennae each attached to a separate receiver branch. The data is demodulated independently for each receiver branch to produce an output signal. The CSI of each output signal is determined and the output signal having the best CSI is chosen for symbol decoding. The decision for choosing which one of the receivers output signals to decode may be performed on a symbol-by-symbol basis.

PLL Adaptation

In a coherent communication system it is customary to use a circuit which tracks and adjusts to carrier frequency changes, commonly known as phase-locked loops (PLLs). The PLL makes adjustments based upon perceived error. If the channel quality is poor, incorrect updates are fed to the PLL and it can become unstable. This can be prevented by using CSI information to adjust the PLL adaptation according to channel conditions. A similar argument holds for phase adjustment loops (PALs) employed in maximum likelihood sequence estimation (MLSE) decoders which track the carrier phase.

Errors and Erasure Decoding

Errors and Erasure decoding is typically employed with forward error correction schemes. Each decoded bit or symbol is identified as good or corrupted. The corrupted bit/symbols are marked as such. A good bit/symbol is detected as is, but an corrupted bit/symbol is ignored. The redundant error correction scheme at a receiver derives a replacement bit/symbol for each corrupted bit/symbol from the remaining bits/symbols. If the corrupted bit/symbol was employed in decoding, then a mistake could have been made which in turn would cause further errors. Hence, erasure decoding enhances bit error rate (BER) performance considerably. To effect such erasure decoding a symbol by symbol estimate of CSI is required.

Currently there is a need for a suitable measure of channel quality (CSI) from the phase angle of the received signal.

SUMMARY OF INVENTION

Signals transmitted through a channel are received, and the channel quality, referred to as a signal-to-impairment ratio (SIR) is determined for the channel. This SIR takes into account not only added white Gaussian noise, but data noise introduced by intersymbol interference, adjacent channel interference, and co-channel interference. An instantaneous and averaged SIR is determined from received signal phase angles, using the principles of ML (maximum likelihood) and MAP (maximum a posteriori) estimation. The SIR values are used in a number of applications. Specifically, they are useful in errors and erasure decoding, post-detection antenna diversity, and coherent detection, employing phase-locked loops.

OBJECTS OF THE INVENTION

It is an object of the present invention to provide a metric for estimating channel state information (CSI) that can be derived from received signal phase angles.

Another object of the present invention is to provide a method of determining CSI that does not require sounding or pilot signals.

Another object of the present invention is to provide a method of determining CSI in systems where only symbol phase angles are available such as systems employing hard-limited signals or systems employing phase detectors.

BRIEF DESCRIPTION OF THE DRAWINGS

The features of the invention believed to be novel are set forth with particularity in the appended claims. The invention itself, however, both as to organization and method of operation, together with further objects and advantages thereof, may best be understood by reference to the following description taken in conjunction with the accompanying drawings in which:

FIG. 1 is a simplified block diagram of a Post Detection Selection Diversity demodulator having two demodulation branches employing signal-to-impairment ratio estimation according to the present invention.

FIG. 2 is a simplified block diagram of a second order phase-locked loop with adaption control employed in a demodulator employing signal-to-impairment ratio estimation according to the present invention.

FIG. 3 is a simplified block diagram of a first embodiment of a demodulator employing signal-to-impairment ratio estimation for errors and erasures decoding, according to the present invention.

FIG. 4 is a simplified block diagram of a second embodiment of a demodulator employing signal-to-impairment ratio estimation for errors and erasures decoding, according to the present invention.

DETAILED DESCRIPTION OF THE INVENTION CSI Metric Derivation

Digital information desired to be transmitted, message data, is encoded and transmitted through a channel to a receiver and heterodyned to a baseband signal. The received baseband signal is sampled into a stream of received samples r_(k), which can be represented as:

    r.sub.k =g.sub.k e.sup.jθ.sbsp.k +I.sub.k e.sup.jΨ.sbsp.k +η.sub.k                                              (1)

corresponding to the k^(th) sampling instant, where g_(k) is the amplitude of the desired signal, and θ_(k) is the corresponding phase angle, I_(k) is the composite amplitude of the interfering signals, Ψ_(k) the corresponding phase angle of the interfering signal, and η_(k) is the additive white Gaussian noise (AWGN) component. It is assumed that one sample per symbol interval is sufficient for estimating the CSI.

From Equation (1) it is observed that the received signal comprises the desired signal and two sources of interference, namely, data noise from interfering signals and AWGN. The data noise may arise when a previously transmitted symbol interferes with a presently transmitted symbol, known as intersymbol interference (ISI). Signals from transmitters geographically adjacent to the receiver, and transmitted on frequencies spectrally adjacent those received by the receiver may be a second source of data noise known as adjacent channel interference (ACI). A third source of data noise is from a transmitter transmitting on the same frequency, but from a communicating cell which is not adjacent to the present cell. This is known as co-channel interference (CCI).

Data noise can be assumed to be a zero mean complex Gaussian random variable similar to the AWGN. This is a practical approximation which simplifies the derivation of the CSI estimate from the received phase angle. Let the variance of the data noise and AWGN be denoted as 2σ_(D) ² and 2σ_(N) ² respectively.

The measure of CSI is defined to be the ratio of the instantaneous signal power g_(k) ², to the variance of the data noise plus the variance of the AWGN component as follows: ##EQU1## where γ_(k) is the instantaneous signal-to-impairment ratio (SIR). SIR is a better indicator of channel quality than signal-to-noise ratio since it also takes into account the influence of data noise which is a function of the channel. If there is no data noise (no ISI, ACI or CCI), then Equation (2) reduces to an instantaneous signal-to-noise ratio.

The problem that remains is to determine γ_(k) which will be used as an estimate of CSI. From Equation (1), the received signal phase angle is: ##EQU2## where η_(ik) and η_(qk) are the in-phase and quadrature components of the AWGN, namely, η_(k) =η_(ik) +jη_(qk).

If the data noise and AWGN are grouped together as a composite complex Gaussian noise component (ξ_(qk), ξ_(ik)) with zero mean and a variance of 2(σ_(D) ² +σ_(N) ²), it may be expressed as: ##EQU3## where ξ_(qk) =I_(k) sin Ψ_(k) +η_(qk) and ξ_(ik) =I_(k) cos Ψ_(k) +η_(ik), are zero mean Gaussian random variables with variance σ_(D) ² +σ_(N) ². A conditional probability density function being the probability of receiving a phase angle (φ_(k) given θ_(k) and γ_(k) is described by: ##EQU4## where erfc is a complementary error function defined by: ##EQU5##

For γ_(k) >>1 and ##EQU6##

Eqs. (5a) and (5b) may be approximated as: ##EQU7##

The above approximation is valid for reasonable SIRs (>3 dB) and for ##EQU8##

An estimate of γ_(k) or CSI may now be obtained employing the maximum likelihood (ML) estimation rule which is expressed as a solution to the equation: ##EQU9##

ln denotes the natural logarithm. Substituting Eq. (6) in Eq. (7) and solving yields: ##EQU10## where γ_(k) ^(ML) is an ML estimate of the SIR, γ_(k).

In many instances the estimate provided for by Eq. (8) will be noisy and may be improved upon by averaging over a number of instantaneous estimates. This modification is useful only if the data samples are independent. Assuming the noise samples are independent from sample to sample and that the channel is quasi-stationary, then N samples may be used to obtain a better ML estimate as a solution to: ##EQU11##

In Eq. (9) it was assumed that samples are required at the same point in each symbol interval and that the channel is stationary therefore γ_(k) is independent of the index k and γ_(k) =γ. Now solving Eq. (9) using Eq. (6): ##EQU12##

The estimate given by Eq. (10) holds only if the channel is quasi-stationary and if there is independence of samples. This is typically true for an indoor portable radio channel.

For a Rayleigh fading channel, an a posteriori estimate may be made for the probability density function for γ_(k). An alternate estimate, being a "maximum a priori" (MAP) estimate may be used to replace the ML estimate, but only for Rayleigh fading channels. The MAP estimate γ^(MAP) is given by a solution to: ##EQU13## with the probability density function P(γ_(k)) given by: ##EQU14## where γ is the average signal-to-noise ratio (SNR) which may be calculated by conventional methods. Now, solving Eq. (11) employing Eqs. (6) and (12) results in: ##EQU15## at the kth sampling point.

Once again, for the quasi-static channel with independent samples, a time averaged estimate may be obtained as a solution to: ##EQU16##

Solving Eq. (14) with Eqs. (6) and (12): ##EQU17##

Equations (8), (10), (13) and (15) provide 4 metrics which may be used to estimate the SIR which can be used as an indicator of channel quality. These may be compared to a minimum threshold value or the inverse of the SIR values may be compared to a maximum value to test for signal quality. The latter comparison does not require a division which the former would, resulting in fewer calculations.

A smoothed metric may be defined as follows: ##EQU18## where λ is a value between 0 and 1 and is typically chosen as a compromise between smoothing the noisy estimate and tracking fast variations of channel quality or SIR. This may be computed recursively as follows: ##EQU19## and ρ_(k) is given by:

    ρ.sub.k =λρ.sub.k-1 +1                      (18)

and ρ₀ =0.

Another, useful metric is to chose γ as the median value of a set of estimates given by:

    γ.sup.med =median{γ.sub.1,γ.sub.2, . . . ,γ.sub.N }(19)

where γ_(k), k=1,2, . . . , N may be computed from Eqs. (8) or (13).

While the use of γ_(k) (the instantaneous SIR) for post detection diversity combining yields improvement in a fading channel compared to differential detection alone, the estimate of instantaneous SIR will be noisy. The performance of post detection selection diversity can be improved by filtering the selection metric with metric filters, if the noisy symbols are independent from symbol to symbol and the channel is quasi-stationary over a duration of a metric filter impulse response. Let γ_(k) be the instantaneous SIR metric at time kT, where T is the symbol interval and k is an integer. Then, a simple filtered selection metric can be written as: ##EQU20## where h(.) are coefficients of the impulse response of a metric filter and L is the length of the impulse response in symbol intervals, chosen to be the average fade duration of the channel.

Eq. (20) represents causal filtering of selection metric γ_(k). If signal samples and instantaneous SIR estimates are stored in memory, the selection metrics can be smoothed by non-causal filtering. This non-causally filtered metric can be expressed as: ##EQU21## where ##EQU22## and x is a floor function representing a largest integer less than x; and x is a ceiling function representing a smallest integer greater than x.

Summary of CSI Metrics

Eight metrics have been derived for CSI based on a differential phase between the received signal phase φ_(k) and the transmitted phase angle θ_(k). These metrics are as follows: ##EQU23## where γ_(k) may be computed from Eqs. (8) or (13). The metric described by Eq. (26) can be computed recursively for a desired number of repetitions to result in a smoothed SIR estimation according to: ##EQU24## and ρ_(k) is given by:

    ρ.sub.k =λρ.sub.k-1 +1                      (27b)

and ρ₀ =0, and λ is a value between 0 and 1.

Metric 6: Median Estimate

    γ.sup.med =medial{γ.sub.1,γ.sub.2, . . . ,γ.sub.N }(28)

where γ_(k), k=1,2, . . . ,N may be computed from Eqs. (8) or (13). ##EQU25## where h(i) are coefficients of the impulse response of a metric filter and L is the length of the impulse response in symbol intervals, chosen to be the average fade duration of the channel. ##EQU26## and x is a floor function representing a largest integer less than x; and x is a ceiling function representing a smallest integer greater than x.

Application of CSI Metrics

The CSI metrics described in Eqs. (22)-(30) may be employed to enhance the performance of a digital radio receiver in many applications such as the following.

a. Post-Detection Selection Diversity

FIG. 1 illustrates a schematic of a digital receiver having 2 branches, "a" and "b", employing post detection selection diversity. Antennae 1a and 1b sense a radio signal which has digital information encoded in it. Front end units 2a, 2b convert the signal sensed by antennae 1a, 1b, respectively, to an intermediate frequency signal and then to a baseband signal which is sampled to produce the stream of samples r_(k) ¹,r_(k) ². Phase detectors 3a, 3b detect a plurality of phase angles φ_(k) ¹,φ_(k) ² from the samples of front end units 2a, 2b, respectively.

SIR calculation units 4a, 4b determine SIRs γ¹,γ² using phase angles φ_(k) ¹,φ_(k) ² and the symbol constellation phase points {θ} stored in a suitable storage means, such as in a read-only memory (ROM) 4c of FIG. 1, according to Eqs. (22)-(30).

A selection logic circuit 5 compares the SIR metrics γ¹,γ², provided by SIR calculation units 4a, 4b, and causes a branch selection device 6 to connect the receiver branch corresponding to the SIR value indicating greater channel quality to a decoder 7. Decoder 7 then decodes the phase angle φ_(k) ^(j) into decoded message bits. This selection is done on a symbol-by-symbol basis. Instead of using a maximum value of γ^(j) decision rule, selection logic unit 5 may employ a minimum value of 1/γ^(j) for receive branch selection decisions. The latter rule avoids a division operation, as indicated by Eqs. (22)-(30).

SIR calculation units 4a, 4b may also employ metric filter units which filter the SIR metrics according to equations (29) and (30). Eqs. (22), (24), (26), (29) and (30) describe metrics best suited for a mobile radio channel, while Eqs. (23), (25), (26) and (28) describe metrics best suited for indoor portable radio channels.

If Eqs. (22) or (24) are used as the metric of channel quality, the metric calculation may be simplified as follows. Let γ¹,γ² represent the SIRs estimated for the output signals from phase detectors 3a, 3b, respectively.

The comparison of γ_(k) ¹ to γ_(k) ² is equivalent to the comparison of |θ_(k) ¹ -φ_(k) ¹ | to |θ_(k) ² -φ_(k) ² | where θ_(k) ¹,φ_(k) ¹ refer to a transmitted symbol phase, and a received symbol phase, respectively, from phase detector 3a; and θ_(k) ²,φ_(k) ² refer to a transmitted symbol phase, and a received symbol phase, respectively, from phase detector 3b.

b. PLL Adaptation

Conventional communication systems typically employ PLLs as shown in FIG. 2 for automatic frequency control (AFC) and for a coherent demodulation. A PLL is an estimator of received signal phase or frequency and employs a closed loop adaptation procedure to track the signal phase or frequency. An antenna 6 receives a radio signal having digital information encoded in it. A front end 8 heterodynes the received carrier signal to a baseband signal which is then sampled. A phase detector 10 extracts phase angles φ_(k) from the baseband signal samples. A set of predetermined phase angles θ is stored in a storage device such as read only memory (ROM) 12. A summer 14 combines received signal phase angles φ_(k) with a corresponding predetermined phase angles θ_(k) and a feedback signal to produce sent to gain units 20a, 20b. Summers 18a, 18b and delay units 16a, 16b are connected as shown in FIG. 2 to utilize the outputs of gain units 20a, 20b to create to the feedback signal and comprise a PLL. The PLL produces a frequency (phase) adjustment ΔΦ which is used to adjust the frequency being received by a receiver closer to the transmitter frequency. An SIR calculation unit 22 estimates SIR values γ_(k) according to Eqs. (22)-(30). For this application, instantaneous estimates of SIR are too noisy to be reliable and hence the metrics given by Eqs. (23), (25), (26) and (28) are recommended. Under poor channel conditions, the PLL becomes unstable, resulting in what is termed loop "hang-up". This can be mitigated by not allowing the PLL to adapt during poor channel conditions, determined by CIR estimates of Eqs. (22)-(30). During the poor channel condition periods, the PLL does not adjust its current carrier frequency estimation and "coasts". A loop adaption control unit 24 monitors the SIR γ_(k). If the value of γ_(k) falls below a predetermined threshold level, the loop switches S1 and S2 are opened and the loop will coast providing a frequency correction corresponding to the frequency estimate obtained prior to the opening of the switches S1 and S2. When the SIR γ_(k) rises above the predetermined threshold value, S1 and S2 are closed allowing for the PLL to function normally.

c. MLSE Phase Decoder

Since initially known sync signals were transmitted as phase angles, they may be compared with corresponding received symbols during their transmission and the metrics Eqs. (22)-(30) can be used as described previously. In the absence of such explicit knowledge of the transmitted phase angles θ_(k), from a set of symbol constellation phase points {θ}, decision directed estimation will have to be used. That is, a tentative decision on θ_(k) is used to estimate the SIR γ_(k). The tentative decision on θ_(k) is provided using a "maximum likelihood" criterion, based on |θ_(k) -φ_(k) |. In the case of a maximum likelihood sequence estimator (MLSE) type decoder, this estimation is performed by estimating the maximum likelihood path of a trellis, along which the value of θ_(k) is most reliable. For a symbol by symbol detecter, this estimation corresponds to using the value of |θ_(k) -φ_(k) | that is the smallest for all possible transmitted set of symbol constellation phase points {θ}.

A MLSE phase decoder for QPSK signals may be implemented using a 4 state Viterbi algorithm using Eq. (5) as the branch metric. From Eq. (5) it is observed that the SIR, γ_(k), is required to be known to compute the branch metric. If a fixed value of γ_(k) is employed, the implementation would be sub-optimal, and this realization results in poor BER. Hence, it is important to use an estimate of γ_(k), in the branch metric calculation. This is equivalent to having an adaptable branch metric that varies with the channel quality.

The MLSE phase decoder is a coherent decoder and requires tracking of the signal phase. This is done through a PLL. Hence, as mentioned above, the PLL adaptation will also be controlled through the estimate of γ_(k).

The estimates provided by Eqs. (22), (26), (29) and (30) are recommended for use with the MLSE phase decoder for branch metric calculation and PLL adaptation in a mobile or indoor portable radio channel.

d. Errors and Erasure Decoding

FIG. 3 illustrates a simplified block diagram of a system employing errors and erasure decoding. In this system, data desired to be transmitted, message data, is encoded with error correction code to add redundancy. The encoded data is mapped to a series of symbols, θ_(k), each chosen from a set of signal constellation points {θ}. These symbols, θ_(k), are transmitted through the channel. The received signal is processed by a front end unit 27 which heterodynes and samples the received signal into a stream of samples, r_(k). A phase detector 28 receives samples r_(k) and determines phase angles φ_(k) from the samples, which are stored in a buffer 31.

An error correction decoder 29 reads the phase angles φ_(k) from buffer 31 and decodes them into a block of symbols or bits ("decoded message data"), e_(k). The decoded message data, e_(k), may be held in correction decoder 29 or placed at a different location in buffer 31.

A signal-to-impairment ratio (SIR) calculation unit 32 determines an SIR value, γ_(k), from the phase angles φ_(k), either from phase detector 28 or from buffer 31. The SIR values correspond to the message data, e_(k), from the phase angles φ_(k) ¹ k

An error detection unit 35 receives the message data, e_(k), and determines if the decoded data block is acceptable. There are many different methods of determining if the data block is acceptable such as calculating parity bits, or matching the decoded block to a look-up table of valid data blocks by computing a difference metric. If the block of decoded message data is acceptable, the block is provided to an output device 37 to be further utilized.

However, if the block is unacceptable, the error detection unit 35 causes an erasure control unit to determine a phase angle in the buffer which has not yet been marked having the lowest SIR value, γ_(k), and marking it as corrupted. Error correction decoder 29 decodes the phase angles φ_(k) which have not been marked, ignoring the marked phase angles, into another decoded block of message data. Correction decoder 29 takes advantage of the redundancy built into the encoding scheme to process the block of message data with an error correction scheme to replace missing data pertaining to marked phase angles φ_(k). The data block is again checked by error detection unit 35 to determine if it is acceptable. If it is, the block of message data is provided to output device 37; if it is unacceptable, the phase angle corresponding to a lowest SIR value which has not yet been marked, is marked as corrupted by erasure control unit 33, and the processing is repeated with all phase angles which have been marked up until this point being omitted from decoding. A limit is set for the number of symbols which may be omitted from the processing. This limit is determined from the error-correcting capability of the code. If the data block is not acceptable when this number of have been omitted, alternate strategy may be employed.

In FIG. 4, a second embodiment is shown. Elements 27, 28 and 32 operate as described in the first embodiment. The SIR calculation unit 32, receives the phase angles φ_(k) and determines a SIR φ_(k) corresponding to each phase angles φ_(k) and passes the SIR to the erasure control unit 33. The erasure control unit 33 marks phase angles φ_(k) in buffer 31 pertaining to phase angles φ_(k) having an SIR estimate below a predetermined threshold value. Error correction decoder 29 takes advantage of the redundancy built into the encoding scheme to process the symbol phase angles φ_(k) with an error correction scheme to replace marked symbols, and reverse map the resulting symbols into message data which is provided to output device 37.

The metrics provided by Eqs. (22) and (24) may be used in either embodiment in determining which symbols to mark as corrupted for erasure with the threshold being selected from Eqs. (26) and (28). Erasure decoding may be implemented in channels which do not employ radio transmission by eliminating the so antenna and front end 27, and connecting directly to a data line.

e. Sample Timing Selection

Another useful application of the metrics given by Eqs. (22)-(30) is sample timing selection. Each symbol of a sampled received signal comprises a number of samples. The object is to select a sample index that maximizes the SIR and decode the samples for each symbol pertaining to the chosen index. The metrics of Eqs. (22)-(30) may be used to select an effective sample within each symbol for the purpose of decoding with minimum error. An averaging type of metric such as Eqs. (23), (25), (26), (29) and (30) accumulated over several symbols is suggested for this application.

The present invention requires knowledge of the transmitted signal phase θ_(k). This is obtained by employing one of the following schemes. The metrics of Eqs. (22)-(30) estimate the SIR in terms of the difference in absolute phase angles between the transmitted signal (θ_(k)) and the received signal (φ_(k)). This is due to Eqs. (5a) and (5b) being expressed in terms of this absolute phase difference is Eqs. (5a) and (5b) may be written in terms of differential phase angles, namely, Δθ_(k) =θ_(k) -θ_(k-1) and Δφ_(k) =φ_(k) -φ_(k-1) corresponding to the transmitted and received signal differential phase angles, respectively. This is an approximation and is obtained by neglecting the cross product of noise terms between the two signaling intervals. With this approximation, the standard ML and MAP estimations as described above and the metrics of Eqs. (22)-(30) may be written in terms of Δθ_(k) -Δφ.sub. k, instead of θ_(k) -φ_(k) This modification results in a SIR that is 3 dB worse than if the estimate was made in terms of absolute phase angle differences, since noise from two symbol intervals is included in computing the differential phase angles. However, since the metrics of Eqs. (22)-(30) are used only in a relative sense (comparison against one another), the use of differential angles for SIR estimation is still effective. The modified metric, using differential phase angles, can be directly employed with differential detectors.

The present invention may be employed in mobile and portable radio receivers that will be used in indoor and mobile radio channels to effect post-detection selection diversity with differential detection, erasure decoding, PLL adaptation (AFC), MLSE phase equalization and sample timing estimation. The technique, however, is general and may be employed in radio receivers for other channels and other applications where an estimate of CSI can be used to enhance receiver performance.

While several presently preferred embodiments of the present invention have been described in detail herein, many modifications and variations will now become apparent to those skilled in the art. It is, therefore, to be understood that the appended claims are intended to cover all such modifications and variations as fall within the true spirit of the invention. 

What is claimed is:
 1. A selection combining receiver for decoding message bits transmitted as symbol phase angles each selected from a set of symbol constellation phase points {θ} in a signal through a channel comprising:a) a plurality of receive branches, each branch comprising:i. an antenna "i" for receiving the transmitted signal; ii. a front end for converting the received signal to an intermediate frequency signal and converting the intermediate frequency signal into a stream of baseband signal samples; iii. a phase detector for detecting a plurality of phase angles φ_(k) ^(i) from the samples; iv. a signal-to-impairment ratio (SIR) calculation unit for determining an SIR value from the detected phase angles φ_(k) ^(i) described by the equation: ##EQU27## where θ_(k) is a transmitted phase angle from the set of symbol constellation phase points {θ} which is closest to φ_(k) ^(i) and N is a number of symbols over which the SIR is desired to be averaged; b) a decoder adapted for decoding the detected phase angles φ_(k) ^(i) into message bits; c) a branch selection device capable of connecting a receive branch to the decoder; d) a selection logic unit coupled to the SIR calculation units and the branch selection device for determining from the SIR values from each receive branch, the receive branch pertaining to the best channel quality and causing the branch selection device to connect the decoder to this receive branch causing the signal from that branch to be decoded.
 2. A method of determining channel quality of a carrier modulated signal encoding digital information as symbol phase angles each selected from a set of symbol constellation phase points {θ} transmitted through a AWGN channel comprising the steps of:a) receiving, heterodyning, and sampling the carrier modulated signal to produce a plurality of baseband received samples r_(k) ; b) decoding the received samples r_(k) into a plurality of phase angles φ_(k) ; c) computing SIR value γ_(k) from phase angles φ_(k) according to the equation: ##EQU28## where θ_(k) is a transmitted phase angle from the set of symbol constellation phase points {θ} which is closest to φ_(k) ; and d) computing a channel quality from the SIR values.
 3. The method of determining channel quality of a carrier modulated signal of claim 2 wherein the equation by which the SIR values are computed is replaced by the equation; ##EQU29## where γ^(i) is the SIR value and is independent of the index k, θ_(k) is a transmitted phase angle from the set of symbol constellation phase points {θ} which is closest to φ_(k), and N is a number of symbols over which the SIR values are desired to be averaged.
 4. The method of determining channel quality of a carrier modulated signal of claim 2 wherein the equation by which the SIR values are computed is replaced by the equation: ##EQU30## where θ_(k) is a transmitted phase angle from the set of signal constellation phase points {θ} which is closest to φ_(k) and γ is an average signal-to-noise ration.
 5. The method of determining channel quality of a carrier modulated signal of claim 2 wherein the equation by which the SIR values are computed is replaced by the equation: ##EQU31## where γ^(i) is the SIR value and is independent of the index k, θ_(k) is a transmitted phase angle from the set of signal constellation phase points {θ} which is closest to φ_(k), is an average signal-to-noise ratio and N is a number of symbols over which the SIR values are desired to be averaged.
 6. The method of determining channel quality of a carrier modulated signal of claim 2 or claim 4 further comprising the step of: exponentially smoothing the SIR values, γ_(k) ^(i), to result in smoothed SIR estimates, γ_(k) ^(i)(s), according to: ##EQU32## where θ_(k) is a transmitted phase angle from the set of signal constellation phase points {θ} which is closest to φ_(k), λ is a value between 0 and 1 with values closer to 1 providing greater smoothing at the expense of faster tracking, and values closer to 0 providing faster tracking at the expense of lesser smoothing.
 7. The method of determining channel quality of a carrier modulated signal of claim 6 which can be computed recursively as follows:a) setting a smoothed signal-to-impairment ratio (SIR) for k=0, γ₀.sup.(s), to a predetermined initial value; b) setting a smoothed smoothing factor ρ₀ to 0 for k=0; c) setting k to 1, and λ to a predetermined value from 0 to 1; d) computing ##EQU33## e) computing ρ_(k) given by: ρ_(k) =λρ_(k-1) +1 f) incrementing k; g) repeating steps "d" through "f" for k less than a predetermined maximum value to reset in smoothed SIR values γ_(k).sup.(s).
 8. The method of determining channel quality of a carrier modulated signal of claim 2 or claim 4 further comprising the steps of:a) calculating the SIR values , γ_(k) ^(i), for K=1,2, . . . ,N to produce {γ₁,γ₂, . . . ,γ_(N) }; b) determining a median value of the calculated SIR values {γ₁,γ₂, . . . ,γ_(N) } being a median SIR γ_(k) ^(med) ; c) employing the median SIR γ_(k) ^(med) as the SIR values.
 9. The method of determining channel quality of a carrier modulated signal of claim 2 further comprising, before the step of computing a channel quality, the step of:filtering the SIR values γ_(k) ^(i) with a metric filter according to the equation: ##EQU34## where h(.) are coefficients of an impulse response of the metric filter and L is the length of the impulse response of the metric filter in symbol intervals, chosen to be the average fade duration of the channel, and employing the filtered SIR values γ_(k)(ave)^(i) in subsequent steps in place of SIR values γ_(k) ^(i).
 10. The method of determining channel quality of a carrier modulated signal of claim 2 further comprising, before the step of computing a channel quality, the step of:filtering the SIR values γ_(k) ^(i) with a metric filter according to the equation: ##EQU35## where and x is a floor function representing a largest integer less than x, x is a ceiling function representing a smallest integer greater than x, h(i) are coefficients of an impulse response of the metric filter and L is the length of the impulse response of the metric filter in symbol-intervals, chosen to be the average fade duration of the channel, and employing the filtered SIR values γ_(k)(ave')^(i) in subsequent steps in place of SIR values γ_(j) ^(i).
 11. The method of determining channel quality of a carrier modulated signal of claim 2 wherein the channel quality is computed by comparing the SIR values to a minimum threshold value with an SIR value above the threshold value indicating an acceptable channel quality, and an SIR value not above the threshold value indicating a channel quality that is not acceptable.
 12. The method of determining channel quality of a carrier modulated signal of claim 2 wherein the channel quality is computed by comparing the inverse of the SIR values to a maximum threshold value, with an inverse SIR value below the threshold value indicating an acceptable channel quality, and the SIR value not below the threshold value indicating a channel quality that is not acceptable.
 13. The method of determining channel quality of a carrier modulated signal of claim 2, 3, 4 or 5 wherein phase angles φ_(k), θ_(k), are differential phase angles defined by:φ_(k) =φ'_(k) -φ'_(k-1), θ_(k) =θ'_(k) -θ'_(k) -θ'_(k-1), where φ'_(k), θ'_(k) are absolute phase angles transmitted at time `k`, and φ'_(k-1), θ'_(k-1) are absolute phase angles transmitted at time `k-1` which precedes time `k`, respectively.
 14. A method of decoding digital information redundantly encoded as symbols from a signal constellation in a signal transmitted through a channel comprising the steps of:a) receiving, heterodyning, and sampling the transmitted signal to produce a plurality of baseband received samples r_(k), with `k` being an index applying to a kth sampling instant; b) determining phase angles φ_(k) ; c) identifying each phase angle φ_(k) as not marked; d) calculating signal-to-impairment ratio (SIR) values from the determined phase angles φ_(k) indicating channel quality according to the equation: ##EQU36## where θ_(k) is a transmitted phase angle from the set of symbol constellation phase points {θ} which is closest to φ_(k) ; and e) computing a channel quality from the SIR values; f) decoding phase angles φ_(k) which have not been marked into a block of message data employing an error correction scheme which takes advantage of the redundancy in the encoding; g) determining if data block is acceptable, andi. if it is acceptable, providing the message data block to a digital output device, and ii. if it is not accentable, identifying a phase angle φ_(k) pertaining to a lowest SIR value which has not yet been marked as marked; and h) repeating steps "f"-"g" until the data block is acceptable for a predetermined maximum number of repetitions.
 15. A method of decoding digital information redundantly encoded as symbols from a signal constellation in a signal transmitted through a channel comprising the steps of:a) receiving, heterodyning, and sampling the transmitted signal to produce a plurality of baseband received: samples r_(k) ; b) determining phase angles φ_(k) from the samples r_(k), with `k` being an index applying to a kth sampling instant; c) estimating which symbol from the signal constellation transmitted corresponds to phase angles φ_(k) ; d) calculating signal-to-impairment ratio (SIR) values for estimated phase angles φ_(k) indicating channel quality according to the equation: ##EQU37## where θ_(k) is a transmitted phase angle from the set of symbol constellation phase points {θ} which is closest to φ_(k) ; and e) computing a channel quality from the SIR values; f) determining for each of the SIR values, if the SIR value exceeds a predetermined value, then erasing the symbol corresponding to the SIR value; g) decoding each symbol which has not been erased by reverse mapping each phase angle φ_(k) to a bit string according to the signal constellation employed in encoding the digital information; h) employing an error correction scheme which takes advantage of the redundancy in the encoding and provides bit strings originally encoded in the phase angle φ_(k) that have been marked, from the remaining bit strings; and i) repeating steps "a"-"h" for subsequent signal transmissions.
 16. The method of decoding digital information of claim 14 or 15 wherein the equation by which the SIR values are computed is replaced by the equation: ##EQU38## where γ^(i) is the SIR value and is independent of the index k, θ_(k) is a transmitted phase angle from the set of symbol constellation phase points {θ} which is closest to φ_(k), and N is a number of symbols over which the SIR values are desired to be averaged.
 17. The method of decoding digital information of claim 14 or 15 wherein the equation by which the SIR values are computed is replaced by the equation: ##EQU39## where θ_(k) is a transmitted phase angle from the set of signal constellation phase points {θ} which is closest to φ_(k) and γ is an average signal-to-noise ratio.
 18. The method of decoding digital information of claim 14 or 15 wherein the equation by which the SIR values are computed is replaced by the equation: ##EQU40## where γ^(i) is the SIR value and is independent of the index k, θ_(k) is a transmitted phase angle from the set of signal constellation phase points {θ} which is closest to φ_(k), γ is an average signal-to-noise ratio and N is a number of symbols over which the SIR values are desired to be averaged.
 19. A system for decoding digital information redundantly encoded as symbols from a set of signal constellation points {θ} in a signal transmitted through a channel comprising:a) a front end capable of receiving, heterodyning, and sampling the transmitted signal to produce a plurality of baseband received samples r_(k), with `k` being an index applying to a kth sampling instant; b) a buffer capable of retaining a plurality of phase angles; c) a phase detector coupled to the front end adapted for providing a plurality of phase angles φ_(k) from the baseband received samples r_(k) and storing the phase angles φ_(k) in the buffer; d) an SIR calculation unit capable of calculating a signal-to-impairment ratio (SIR), γ_(k), for each phase angle φ_(k) indicating channel quality according to the equation: ##EQU41## where θ_(k) is a transmitted phase angle from the set of symbol constellation phase points {θ} which is closest to φ_(k) ; and d) computing a channel quality from the SIR value; e) an error correction decoder capable of estimating a block of message dam, e_(k), corresponding to the transmitted signal constellation points {θ} from phase angles φ_(k) which have not been marked as corrupted, and adapted for employing the redundancy inherent in the encoding scheme to calculate message data pertaining to phase angles marked, as corrupted; f) erasure control unit coupled to the SIR calculation unit, for sequentially marking a phase angles φ_(k) having the lowest SIR value of remaining unmarked phase angles as a corrupted phase angle in the, buffer; and g) an output device capable of utilizing the message data, e_(k) ; h) an error detection unit capable of determining if the block of message data, e_(k), is acceptable, andi. if acceptable providing the block of message data bit to the output device, and ii. if the block of message data is not acceptable, then iteratively for a predetermined maximum number of repetitions, causing the erasure control unit to mark a phase angle having the lowest SIR of all phase angles which have not been marked, and causing then causing the error correction decoder to estimate an updated block of message data, e_(k) from the remaining phase angles which have not been marked. 